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Abstract 

We study the two-dimensional extension of the Fano-Anderson model on the 
basis of a two-dimensional optical waveguide array with a bi-modal defect. We 
demonstrate numerically the persistence of the Fano resonance in wavepacket 
scattering process by the defect. An analytical approximation is derived for 
the total scattered light power. 
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Optical waveguide arrays are known to be very efficient experimental realizations of one- 
(1D) and two-dimensional (2D) discrete systems^. Recent achievements both in sample 
preparation techniques and in experimental methods allowed for observation of several im- 
portant phenomena known from linear and nonlinear lattice theory^. These experimental 
successes, in turn, stimulated theorists in considering other possible implementations of 
waveguide arrays, being of potential interest both in fundamental and applied science. One 
of such promising directions of research, being extensively developed during recent years, is 
related to the possibility of direct observation of the Fano resonance- in processes of wave 
scattering by nonlinear localized excitations^ and in the framework of different linear and 
nonlinear generalizations- of the Fano- Anderson model^^. In all such models a defect (ei- 
ther geometrical or dynamical) is locally attached to an otherwise translationary invariant 
ID waveguide array. This defect mediates one or several additional propagation paths for 
waves, locally increasing the dimensionality of the system. The Fano resonance manifests 
through the resonant reflection of light (as the result of destructive interference between the 
propagation paths) at a particular angled, defined by the defect parameters. 

In this letter we consider the two-dimensional extension of the Fano-Anderson model. 
The setup is based on a 2D waveguide array, see Fig. ^a). An additional defect is assumed 
to be locally coupled to the central waveguide {0,0}. Experimentally it can be realized 
e.g. by making the central waveguide to be bi-modal, so that the defect site corresponds to 
the second mode of this waveguide. Assuming that the overlap between the two modes of 
the bi-modal waveguide is much stronger than the overlap between the second mode of the 
central waveguide and the first (principal) mode of an adjacent waveguide, the coupling to 
the defect is essentially local. 

Within the framework of coupled-mode theory^, the effective model equations describing 
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the evolution of the {n, m}th waveguide mode amplitude u n) m and the defect mode (the 
second mode of the central bi-modal waveguide) amplitude ip along the propagation distance 
z can be written as*- 

i + VAu ntm + ei])8 nfi 8 mfl = , 
dib 

i^ + E^j + ew ,o = , (1) 
dz 

where Au nm = (u n+ x im + tt n _i, m + "Un,m+i + w niTO _i), V is the effective coupling between 
modes of the adjacent waveguides, e is the effective coupling between the two modes of 
the central waveguide, and E describes the difference between the propagation constants of 
the two modes. As usual for Schrodinger-type equations, Eqs. (0) conserve the total power 
(norm) P 

P = J2Pn,m+ Pdef = ^2\u„, m \ 2 + \lp\ 2 , (2) 
n,m n,m 

directly related to the electric field power of a light beam propagating in the array. 

Considering stationary solutions u n>m {z) = A n>m exp [iuz] and ipo = B exp [icoz], we obtain 
the following equations for A n>m : 

e 2 

uA nm = VAA nm + pr^4o,0^n,0^m,0- (3) 

(u-E) 

Far away from the defect site {0, 0} eigenstates of the system © asymptotically approach 
plane waves A nm = A exp [i(k x n + k y m)} with the dispersion relation and group velocities 
given by 

to = 2^(008^3; + cos ky), (4) 
v g = —2V(smk x ,smky), (5) 

respectively. The absolute value of the group velocity (0) is shown in the density plot in 
Fig. ^c). It reaches its maximum 



Analogously to the ID Fano-Anderson models, the defect coupled to the central site 
acts as an effective scattering potential for traveling waves with strength depending on the 
spectral parameter (u) of the incoming wave. The scattering potential diverges at u = E, 
this constitutes the resonance condition. Combining this with the dispersion relation 
we obtain the implicit relation 



which defines the resonance condition curve in the parameter space (k x , k y ) plotted in 
Fig. ^c) with thick lines. Thus, for a plane wave coming at a certain angle (i.e. with 
the fixed relation between k x and k y ) one has a unique pair of parameter values (k x , k y ) for 
which the resonance in the scattering process is expected to occur, see Fig. [TJd) . 

We performed numerical simulations of wavepacket scattering in the array of 81 x 81 
waveguides with fixed boundary conditions. The initial wavepacket is taken in the Gaussian 
form A nm (z = 0) = Aexp [i(k x n + k y m)] with a = 0.01, so that the effective half-width of 
the packet in real space is about 15 sites, while the half-width in the reciprocal /c-space is 
Afc « 0.2. The defect parameter E is chosen to be reasonably small, E = 0.1, so that the 
resonance condition curve defined by Eq. (jHJ) in the parameter space (k x , k y ) passes close 
enough to the points \k x \ — \k y \ — ir/2 corresponding to the maximum group velocity of the 
wavepacket. This allows us to save computational time in numerical simulations described 
below. 

In Figs. 12131 characteristic snapshots of the power density p n ^ m © distribution across the 
array, after the scattering process has occurred, are shown for resonant and non-resonant 
cases, respectively. While for the non-resonant case the wavepacket remains practically 
undistorted after passing the defect site, for the resonant case a considerable amount of 
the initial wavepacket power is scattered in different directions. In order to estimate the 




(6) 
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amount of scattered power we define an approximate boundary between "transmitted" and 
"scattered regions" , see dashed line in Figs. 12131 The total scattered power (the power in 
the "scattered region") as a function of the wavepacket number k x for the case of diagonal 
propagation (k x = k y ) is plotted in Fig. 0] One can clearly observe the resonant peak at 
k x fa 1.56, which is in full agreement with the resonance condition k x = arccos[-E/(4\^)] 
obtained from Eq. (JBj). 

The resonant peak power of scattered light, as well as the characteristic resonance width, 
depend on the initial wavepacket parameters and the parameters of the defect. In order 
to have at least a qualitative understanding of the problem, we make a connection to the 
ID Fano-Anderson model by noting that due to the locality of the defect essentially only 
the diagonal row in the wavepacket interact with the defect site, while the rest passes 
through practically undistorted, see Fig. |21 Referring to the exact result for the transmission 
coefficient of plane waves in ID case 



W e 4 + W 2 (u- E) 2 sm 2 q { ) 

(q is the wavenumber of the ID plane wave) and using the Fourier transform of the diagonal 
part of the initial Gaussian wavepacket, we finally arrive to the following expression for the 
scattered power: 

pR(k x ) = J!Lfdq{\l- T^ D \q)] e -(^) 2 /(2«)| ? (8) 
\J iia J L - 1 ' 

where P d is the power of the diagonal row in the initial wavepacket: 

P d =Yl Pn,m(z = 0) « \ff- (9) 
n=m * 

Despite the simplicity of this ID approach, it nevertheless gives us a reasonably good es- 
timate for the scattered power, see solid line in Fig. 0] The noticeable reduction in the 



peak scattered power is due to the fact, that the effective ratio of the original wavepacket 
which essentially interacts with the defect is actually more than one row. Analysing the 
results given by Eqs. (JH1 EJ) 5 w e conclude that the Gaussian parameter a performs a twofold 
influence on the amount of scattered power. First of all, it controls the ratio of the initial 
wavepacket which effectively interacts with the defect, see Eq. ©. The more localized the 
initial wavepacket is in real space, the higher is the ratio of it which interacts with the 
defect. On the other hand, the more it is localized in real space, the less it is localized 
in the reciprocal /c-space. This, in turn, leads to a broadening of the transmission curve 
for pure plane waves, making the effective resonance width larger and, at the same time, 
decreasing the maximum of the scattered power. This broadening effect, however, can be 
controlled by changing the relative strength of the coupling to the defect e/V. Indeed, the 
rate of broadening is defined by the ratio of the wavepacket width in fc-space to the 
effective resonance width Ap, the latter is known for the ID Fano- Anderson mode l 3161 18 to 
be completely determined by the ratio e/V. 

To conclude, we discussed the 2D Fano- Anderson model based on a 2D optical waveguide 
array with a bi-modal central waveguide. Using numerical simulations we demonstrate the 
persistence of the Fano resonance in this model, which should be detectable experimentally. 
On the basis of a simple ID analogue, we derive a reasonably good approximation for the 
total scattered power of light and study its dependence on the initial wavepacket parameters 
and the model parameters. 

R. A. Vicencio's e-mail address is rodrigov@mpipks-dresden.mpg.de 
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List of Figure Captions 

Fig. 1. In a) and b) the scheme of a 2D waveguide array with the defect on the central site, 
and the concept for incoming, reflected, and transmitted waves in the scattering problem 
are shown; c) Density plot of the plane waves group velocity (absolute value), V = 2; 
d) Resonance diagram: Thick solid and dashed lines represent the resonance condition 
curve given by (jUJ) for E = 0.1 and E — 1, respectively. Thin solid and dashed lines 
represent the (k x ,k y ) parameter submanifold corresponding to the uni- directional motion of 
the wavepacket with v y = v x and v y = v x /2, respectively. 

Fig. 2. Power density distribution p mn (J2J) at z — 10 (after the scattering process). 
Initial wavepacket parameters are: k x = k y = 1.55 (resonant scattering), A = 0.1. Coupling 
constants are: V = 2, e = 1. 

Fig. 3. Power density distribution p mn (J2J) at z — 12.5 (after the scattering process). 
Initial wavepacket parameters 2.2 (non-resonant scattering), A = 0.1. 

Coupling constants are: V — 2, e = 1. 

Fig. 4. The total scattered power for the diagonal propagating wavepacket ik x = k y ) 
calculated from numerical simulations (points) and by the approximation © (continuous 
line). All the parameter values are the same as in Figs. 12131 
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Fig. 1. In a) and b) the scheme of a 2D waveguide array with the defect on the central site, 
and the concept for incoming, reflected, and transmitted waves in the scattering problem 
are shown; c) Density plot of the plane waves group velocity (absolute value), V = 2; 
d) Resonance diagram: Thick solid and dashed lines represent the resonance condition 
curve given by (jHJ) for E — 0.1 and E — 1, respectively. Thin solid and dashed lines 
represent the (k x ,k y ) parameter submanifold corresponding to the uni- directional motion of 
the wavepacket with v y = v x and v y = v x /2, respectively. 
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Fig. 2. Power density distribution p mn (J2J) at z = 10 (after the scattering process). Initial 
wavepacket parameters are: k x = k y = 1.55 (resonant scattering), A = 0.1. Coupling 
constants are: V = 2, e = 1. 
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Fig. 3. Power density distribution p mn (J2J) at z — 12.5 (after the scattering process). Initial 
wavepacket parameters are: k x = k y = 2.2 (non-resonant scattering), A = 0.1. Coupling 
constants are: V — 2, e = 1. 
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Fig. 4. The total scattered power for the diagonal propagating wavepacket (k x = k y ) calcu- 
lated from numerical simulations (points) and by the approximation (jSJ) (continuous line). 
All the parameter values are the same as in Figs. 12131 
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